Realizing bright matter- wave soliton collisions with controlled relative phase 
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We propose a method to split the ground state of an attractively interacting atomic Bose-Einstein condensate 
into two bright solitary waves with controlled relative phase and velocity. We analyze the stability of these 
waves against their subsequent re-collisions at the center of a cylindrically symmetric, prolate harmonic trap 
as a function of relative phase, velocity, and trap anisotropy. We show that the collisional stability is strongly 
dependent on relative phase at low velocity, and we identify previously unobserved oscillations in the collisional 
stability as a function of the trap anisotropy. An experimental implementation of our method would determine 
the validity of the mean field description of bright solitary waves, and could prove an important step towards 
atom interferometry experiments involving bright solitary waves. 
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Bright solitary waves (BSWs) in attractively-interacting 
atomic Bose-Einstein condensates (BECs) are an intriguing 
example of a nonlinear wave phenomenon in a degener- 
ate quantum gas CHS. In the mean-field, Gross-Pitaevskii 
equation (GPE) description, BSWs in a quasi- ID BEC with 
no external (axial) trapping potential correspond exactly to 
bright solitons in the focusing nonlinear Schrodinger equation 
(NLSE) 
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This integrable equation describes a diverse range of physical 
systems in addition to BECs |4|, and its bright soliton solu- 
tions have been extensively studied in the context of nonlinear 
optics EHS). The addition of a harmonic trapping potential 
to Eq. ([T]) leads to non-integrability, but BSWs in harmoni- 
cally trapped, quasi- ID BECs remain, in the GPE description, 
highly soliton-like; they collide elastically over a large param- 
eter regime, and their asymptotic trajectories follow particle 
models applicable to solitons (lOl [TT). Relaxing the quasi- 
1D restriction reduces the soliton character of BSWs further, 
but BSWs in the 3D GPE description retain many soliton-like 
characteristics |12], including an absence of dispersion and 
the existence of a well-defined relative phase between BSWs. 
If actual 3D BSWs possess a well-defined relative phase when 
realized experimentally, one can envisage a BSW interferome- 
ter, akin to current matter- wave interferometers but leveraging 
the small size, coherence, and non-dispersive nature of BSWs 
El . Indeed, BSWs have already been proposed as a metrolog- 
ical tool for the study of atom- surface interactions fT3l 

Experiments to date have produced both individual CD and 
multiple ElO BSWs as remnants from the collapse l|14[[T5) of 
a larger condensate. These BSWs were capable, in the case of 
multiple BSWs, of surviving many mutual re-collisions at the 
trap center (2][3]|. In the quasi- ID regime the observed BSW 
motion matches the GPE description of BSWs with relative 
phase O = n El ITS). In the 3D regime, however, BSWs are 
not universally stable against multiple re-collisions; numeri- 
cal integration of the GPE reveals that slow 3D BSWs retain 
their form for fewer collisions when their relative phase, O, 



is equal to than when O = n [12]. The long lifetimes of 
3D BSWs seen in experiment thus seem to imply that their 
relative phase O = n O [T2]|. Modulational instability and 
the shorter lifetime of colliding 3D BSWs when O = have 
been identified as contributory causes to these apparent anti- 
phase relations in both regimes EH OS] 03. However, recent 
simulations of BSW collisions incorporating quantum noise 
have been interpreted as showing the dynamics and collisional 
stability of BSWs to be phase-independent, with the dynam- 
ics for all relative phases corresponding to the GPE descrip- 
tion for <D = 7i ifTHl . Furthermore, no 3D GPE simulation 
of the collapse process has produced BSW remnants match- 
ing those observed in experiment lfT8l . These considerations 
leave the question open: are experimentally observed atomic 
BSWs well-described by an effective single -particle wavefunc- 
tion, propagated by the GPE? 

In this letter we propose an experiment to answer this ques- 
tion. We describe a method which, in a velocity- and phase- 
controlled way, splits a single BSW in an axisymmetric har- 
monic trap into two outgoing BSWs which repeatedly re- 
collide at the trap center. Using the GPE, we demonstrate that 
such pairs of BSWs can be realized in an atomic BEC, and 
analyze their dynamics and collisional stability. We explore 
the crossover from quasi- ID to fully 3D regimes, examining 
the effects of velocity and relative phase on the number of col- 
lisions for which the BSWs remain soliton-like, Cid- In ad- 
dition to demonstrating the expected Cid phase dependence 
for slow collisions, we show that Cid has a strong oscilla- 
tory dependence on the trap anisotropy, due to resonances be- 
tween the frequency of the BSWs' radial oscillations and the 
frequency with which the two BSWs collide. The experimen- 
tal presence (absence) of the predicted phase and anisotropy 
dependencies would indicate the (in) sufficiency of the GPE 
description of BSWs; either of these outcomes would be an 
important result. In particular, sufficiency of the GPE de- 
scription implies a well-defined relative phase between BSWs, 
paving the way for future atom interferometry experiments 
using BSWs, towards which our proposed BSW generation 
method would represent an important step. 

We begin with the GPE for a BEC of N atoms of mass m 
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and (attractive) s-wave scattering length a s < 0, held within a 
cylindrically symmetric, prolate harmonic trap, 
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where g3D = 4nNa s h 2 /m, the condensate wavefunction ^(r, t) 
is normalized to one, and V(r) = m[(jj 2 x 2 + oj 2 (y 2 + z 2 )]/2, 
where oj x and cj r > oj x are the axial and radial trap frequen- 
cies. With strong radial confinement the system is quasi- 
1D and can be described by Eq. ^ with ^(r, t) — > *¥(x, t), 
V(r) — > mco 2 x 2 /2, and — > £id = 2%cj r \a s \N fT9ll . The re- 
sulting configuration has two key length scales; the harmonic 
length ao = ^h/moj x , and the soliton length bo = n 2 '/mgm. 
Rescaling all lengths to units of bo and all times to units of 
fi 3 /mg 2 D 1 10] produces the dimensionless ID GPE 



diHx, t) 



~2dx 2 



*Kx,t), (3) 



where oj = (bo /ao) 2 is a dimensionless effective trap strength 
l20l . In this letter, we consider the effects of abruptly increas- 
ing the scattering length magnitude in such a BEC from initial 
a° s to a s = a 2 a° s (a > 1 and a s ,a® < 0), with initial condition 



iff(x, t = 0) = iff (x) 
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where if/ a (x) is the BSW ground state of the BEC for scattering 
length a° s . We first consider the quasi- ID limit, where a sta- 
ble ground state if/ a (x) always exists, and then 3D, where the 
stable ground state </^(r) exists only for \a° s \ below the critical 
value for the onset of collapse, \a° s \ fl5\ . 

The ground state if/ a (x) may be made by using a magnetic 
Feshbach resonance to adiabatically change the scattering 
length from being initially repulsive to a negative value, a° s , 
with \a° s \ < \a° s \. The rapid change from a° s to a s = a 2 a° s could 
then exploit the same resonance. The density modulation that 
transforms if/ a (x) into i/fo(x) may be achieved with a second 
internal atomic state in an interference protocol: writing the 
total state of the condensed atoms as i//+(x)\+) + if/-(x)\-), we 
begin with all atoms in internal state |+) [i.e., if/+(x) = i// a (x) 
and iff-(x) = 0]. Applying a resonant n/2 pulse to the inter- 
nal state transition yields if/+(x) = if/-(x) = i/r a (x)/ V2. We 
now imprint equal and opposite momenta on the two internal 
states, giving i//+(x) = exp[±i(Kx + (5>)/2]i// a (x)/ V2, which is 
then transformed into 



if/+(x) = cos [(Kx + 0)/2] iff a (x), 
iff-(x) = i sin [(Kx + 0)/2] fc(x), 



(5) 



by a second n/2 pulse. Using resonant light to rapidly ex- 
pel atoms in state |-) from the trap leaves those in state |+> 
with wavefunction \j/o(x) [Eq. Q], with k = a 2 K, and with 
O determined by the phase accumulated at the center of the 
BSW. Note that the loss of atoms between if/ a (x) and if/o(x) is 
balanced by the change in normalization; N denotes the ini- 
tial atom number. There are many potential implementations 
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FIG. 1. (color online). The structure of a multi- soliton pulse. The 
panels show the soliton amplitudes Ay, velocities V}, and fractions 
associated with the initial condition ifro(x) [Eq. {6}] in the NLSE 
[Eq. {T}], computed using a numerical scattering transform l24l . as a 
function of spatial modulation frequency Panels (a-c) correspond 
to a = 2 and (d-f) to or = 2.2. Relative phases are O = (+), 
n/4 ( ), n/2 (A), 3^/4 (□), n (o). Soliton fraction is the ratio of the 
combined norm of the constituent solitons, 2j 2Aj, to the total norm 
£1 \i//o(x)\ 2 dx (23). In the limit k -> oo, when or = 2, A 7 1/8 
[Zy2A,- ^ J" \Mx)\ 2 dx 1/2], and V,- ^ ±^/8 013. 



of this protocol. Using a two-component GPE, we have sim- 
ulated an implementation that uses 85 Rb atoms in the quasi- 
1D regime with an applied magnetic field gradient to transfer 
momentum: we find this simple prototype to be capable of 
generating initial conditions close to Eq. ^ using current ex- 
perimental technology l2T1l . 

Neglecting the axial trapping (setting oj = 0) the ID GPE 
[Eq. ([3])] reduces to the NLSE [Eq. ([lj], and the ground state 
of the BEC before the change in scattering length, i// a (x), is 
a single, stationary bright soliton [23 ]. After density modula- 
tion 



if/o(x) 



=i sech (2^) cos (£ + !)- (6) 



Solutions of the NLSE for this initial condition [Eq. ([6])] are 
well-known in the context of nonlinear optics (7H9) . The case 
k = has been studied analytically by Satsuma and Yajima 
171 using the inverse scattering transform (1ST) [5 ]: for inte- 
ger a = J, Eq. ^ consists of a bound state, or multi-soliton 
pulse, of / solitons with unequal amplitudes Aj and zero ve- 
locity (Vj = 0). For non-integer a = J + /3, Eq. ^ consists 
of / solitons plus radiation, with the norm of the soliton com- 
ponent given by £j 2A ; - JTJ E3- The modulated case (general 
k) has been considered both analytically and numerically by 
Kodama and Hasegawa [8 ] and Afanasjev and Vysloukh |9). 
Fig. [T] shows how the modulation alters the character of a two 
soliton pulse (a > 2): beyond a certain threshold value of k the 
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FIG. 2. (color online). Generation of BSWs with controlled relative 
phase via the interference protocol, in the quasi- ID limit. Panels 
(a-f) show the evolution of ID GPE with trap frequency oj = 0.02 
[oj = inset in (a-d)] and initial condition if/o(x) for a = 2, O = 
and k = (a), £ = 2 (b), = 4 (c), £ = 6 (d), and O = * and 
k = 4 (e), = 6 (f), computed using a pseudospectral split-step 
method. Particle model |10 | BSW trajectories, for effective masses 
and velocities obtained from the numerical scattering transform of 
i//o(x), are overlaid as lines in (c-f). Panels (e) and (f) reproduce (c) 
and (d) for the case O = tt to show the difference in collision profile. 
The density (color) axes are normalized by c = 0.35 (inset c = 0.25) 
in (a) and c = 0.12 (inset c = 0.07) in (b-f). 



pulse "splits" into two solitons with equal amplitudes, oppo- 
site velocities, and relative phase O, plus a negligible radiation 
component. Crucially, control of the modulation corresponds 
to control over the relative velocity and phase of a pair of gen- 
erated bright solitons. 

A similar correspondence exists in the presence of axial 
trapping (oj > 0). In this regime the ID GPE [Eq. ([3])] has 
no soliton solutions so we study the dynamics of initial con- 
dition ifso(x) numerically, concentrating, for simplicity, on the 
case a = 2 (other cases a > 2 are similar except for a slightly 
altered relationship between k and the resulting soliton speed). 
A pair of equal amplitude BSWs are generated with relative 
phase O and velocities controlled by k [Fig. [2j. The axial 
trap confines the outgoing BSWs and causes subsequent re- 
collisions at the trap center, and the relative phase upon re- 
collision is always identical to the original imposed relative 
phase ifTOl . The BSWs remain highly soliton-like: the den- 
sity profile during BSW collisions is similar to that for bright 



FIG. 3. (color online). Stability of BSW collisions in 3D. Panel (a) 
shows the number of ID-like BSW collisions, Cid, as a function of 
k, O, and A |26l . Effective trap frequency oj = 0.02 and a = 2. 
Also shown is the evolution of the positively displaced BSW posi- 
tion, x s (dashed green line, right vertical axis), and the full width at 
half maximum of the integrated radial density distribution, cr r (solid 
red line, left vertical axis), at the indicated points on the k = 4 (b- 
d) and k = 3.25 planes in (a). The computation leverages the radial 
symmetry of the problem, using a pseudospectral split-step method 
in 2D cylindrical coordinates. 



solitons (6) [Fig. [2jd,f)], the BSW trajectories are well de- 
scribed by a particle model (lOl HU [Fig. [2jc-f)], and the 
BSWs are stable against their mutual collisions, insofar as 
they retain their form for a sufficiently large number of col- 
lisions that atom losses, unaccounted for in the GPE, would 
be the lifetime-limiting factor in an experiment. 

Moving beyond the quasi- ID regime, we can generate pairs 
of 3D BSWs with controlled velocity and relative phase using 
the same method. Dynamics in the radial directions can affect 
the stability of the BSWs, however; in certain cases this dras- 
tically reduces the number of collisions for which they retain 
their form. Parameterizing the quasi- lD-to-3D transition by 
the trap anisotropy X - a) r la) x , we write the 3D GPE as 



d<A(r,0 
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2tt, 
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where we use soliton variables [25], and V(r) = oj z [x z + 
A 2 (y 2 + z 2 )]/2. We again study the dynamics of the BSWs 
numerically, quantifying their stability against collisions in 
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terms of their positions and maximum integrated axial den- 
sities at the point of maximum separation — this being much 
easier to measure, on typical experimental scales, than the ex- 
act density profile during the collision. Fig. |3ja) shows how 
the number of ID-like collisions Cw (taken to be those where 
the positions and maximum integrated densities of the BSWs 
subsequently return to within 75% of their original values) 
depends on velocity, relative phase, and trap anisotropy. We 
term these collisions ID-like because all collisions of quasi- 
1D BSWs satisfy these criteria (Cid — > oo). 

As expected, Fig. [3ja) shows that Cid is strongly depen- 
dent on the relative phase at low velocity, with the BSWs be- 
ing most stable around O = n outside the quasi- ID regime 
l26l . At higher velocity this phase-dependence weakens and 
the quasi- ID regime is reached at lower anisotropy. Fig.|3ja) 
also reveals a previously unobserved feature: Cid shows a 
strong, oscillatory dependence on the anisotropy at all veloc- 
ities. This dependence arises from the BSWs being broken 
up by the transfer of energy to radial oscillations [Fig. [3Jb- 
d)]. These oscillations are started by the abrupt change in 
scattering length, and subsequently amplified by collisions if 
the BSWs collide when their radial width is close to its os- 
cillatory maximum. The frequency of the radial oscillations 
is primarily determined by co r \ this leads to the observed os- 
cillations of Cid as a function of A. The amplifying effect of 
collisions decreases with the BSW velocity, and at low ve- 
locity a phase-dependent amplification of the radial oscilla- 
tions emerges [Fig.[3je-g)], which we attribute to the higher 
densities at the point of collision when O = delivering a 
larger "kick" than when <D = n. However, for intermediate 
phases symmetry-breaking population transfer |12]|27) during 
collisions also contributes to the reduction in Cid [Fig. [3jf)]. 
Within the GPE description, Fig. [3] represents a comprehen- 
sive prediction of the BSW dynamics resulting from our split- 
ting protocol. Experimental observation of the dynamics we 
predict would support the validity of the GPE description of 
BSWs and, in the case of the oscillatory dependence of Cid 
on A, open the possibility of controlling the BSW lifetime di- 
rectly. 

To conclude, we have proposed an experiment that pro- 
duces a pair of BSWs with controlled relative phase and veloc- 
ity in a harmonically trapped atomic BEC, and we have ana- 
lyzed the subsequent collisions of these BSWs using the GPE. 
In the quasi- ID regime the BSWs are highly soliton-like and 
stable against their re-collisions. In the fully 3D regime, we 
confirm that the collisional stability of the BSWs depends on 
their relative phase and velocity, and demonstrate for the first 
time a strong oscillatory dependence on the trap anisotropy. 
The presence, or absence, of these effects in experiments pro- 
vides a direct test of whether experimentally observed atomic 
BSWs can be described in terms of a coherent effective single- 
particle wavefunction, propagated by the GPE. 

We thank N. G. Parker, S. A. Wrathmall, and P. M. Sutcliffe 
for many stimulating discussions, and UK EPSRC (Grant No. 
EP/G056781/1), the Royal Society (SLC), and Durham Uni- 
versity (TPB) for support. 



* t.p.billam@durham.ac.uk 
[1] L. Khaykovich et al, Science 296, 1290 (2002). 
[2] K. E. Strecker et a/., |Nature| 417, 150 (2002). 
[3] S. L. Cornish, S. T. Thompson, and C. E. Wieman, Phys. Rev. 

Lett.| 96, 170401 (2006). 
[4] T. Dauxois and M. Peyrard, Physics of Solitons (Cambridge 

University Press, 2006). 
[5] V. Zakharov and A. Shabat, Zh. Eksp. Teor. Fiz. 61, 118 (1971) 

[Sov. Phys. JETP 34, 62 (1972)] 

[6] J. P. Gordon, [OpETeE] 8, 5 96 (1983). 

[7] J. Satsuma and N. Yajima, [Prog. Theor. Phys. Suppl.| 55, 284 

(1974). 

[8] Y. Kodama and A. Hasegawa, |Opt. Lett.|16, 208 (1991). 
[9] V. V. Afanasjev and V. A. Vysloukh,' !^ Opt. Soc. AmT^l ll, 
2385 (1994). 

[10] A. D. Martin, C. S. Adams, and S. A. Gardiner, Phys. Rev. Lett. 
98, 020402 (2007); Phy s. Rev. A 77, 013620 (2008). 

[11] D. Poletti et a/., |Phys. Rev. Lett. 101, 150403 (2008). 

[12] N. G. Parker et al, J . Phys. B 41, 045303 (2008); Phy sicaD] 
238,1456(2009). 

[13] S. L. Cornish et al, Physica D 238, 1299 (2009). 

[14] E. A. Donley et al, Nature 412, 295 (2001); L. P. Pitaevskii, 
|Phys Lett. A| 221, 14 (1996). 

[15] N. G. Parker et al, J. Phys. B 40, 3127 (2007). 

[16] U. Al Khawaj agfq/., Phys. Rev. Lett. 89, 200404 (2002); K. E. 
Strecker et al, N. J. Phys .|5, 73 (2003). 

[17] L. D. Carr and J. Brand, |Phys. Rev. Lett. 1 92, 040401 (2004). 

[18] B. J. Dabrowska-Wuster, S. Wiister, and M. J. Davis, N. J. Phys. 
11, 053017 (2009). 

[19] In this case, ^(r, t) factorizes into ^(x, t) and the radial har- 
monic ground state, ^(r, t) = ^¥(x, t) ^moj r lnfi exp[-maj r (y 2 + 
z 2 )/2h], the latter of which can be integrated out 1 10]. 

[20] Ensuring unit norm for if/(x, t) requires if/(x, t) - ^fb^(x, t). 

[21] We consider an experiment using the hyperfine ground states 
\F = 2,mf = -2) and |3, -2) of 85 Rb. We assume an instanta- 
neous 7i / 2 pulse on the two-photon rf/microwave transition cou- 
pling |2, -2) and |3, -2), and evolve the resulting state for time 
r with magnetic field B = [B + C(x - x )]B. We assume the 
scattering length of the |3, -2) state and the inter-state scatter- 
ing length is aoo = -443a and model that of the |2, -2) state 
as aoo[l - A/(B + Cx - B F )], where A = 10.7 Gauss and 
B F = 155.0 Gauss l22l . In soliton units, the components ex- 
perience potentials T ± y(x - x ), where T = \gF\\^if\/^B(Bo - 
B F )/4ma) 2 \a s \ 2 N 2 and y = \g F \\m f \ti B Ch/8m 2 u 3 r \a s \ 3 N 3 . After 
time r we assume another instantaneous tt/2 pulse and instan- 
taneous expulsion of the |3, -2) component with resonant light. 
An approximation to Eq. ^ is obtained when, for example 
gj x * 10 x 2tt Hz, oj r * 150 X 2n Hz, N » 800, a s * -20a 
(Bo ~ 166 Gauss), C « 75 Gauss cm -1 and r « 10yus. 

[22] N. R. Claussen et a/., |Phys. RevTA[ 67, 060701(R) (2003). 

[23] Note, for any amplitude A, that A sech (Ax) is a stationary soli- 
ton solution to Eq. ([TJ, with norm 2A. 

[24] G. Boffetta and A. R. Osborne, |rComp. Phys.| 102, 252 (1992). 

[25] Ensuring unit norm for if/(r, t) requires if/(r, t) = b 3 Q /2x ¥(r,t). 
Also, in the quasi- ID limit, the integrated axial density 
ff™Jiff(r,t)\ 2 dydz^\iff(x,t)\ 2 

[26] In the case k = 3 quasi- ID behavior is reached (the phase de- 
pendence of Cid ends) at A « 20, but this has been omitted from 
the plotted range for clarity. 

[27] L. Khaykovich and B. A. Malomed, ]Phys. Rev. A 74, 023607 
(2006). 



